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Abstract. We introduce a max-plus analogue of the Petrov-Galerkin finite 
element method, to solve finite horizon deterministic optimal control problems. 
The method relies on a max-plus variational formulation, and exploits the 
properties of projectors on max-plus scmimodulcs. We obtain a nonlinear 
discretized semigroup, corresponding to a zero-sum two players game. We 
give an error estimate of order V At + Ax(At)~ 1 , for a subclass of problems 
in dimension 1. We compare our method with a max-plus based discretization 
method previously introduced by Fleming and McEneaney. 



1. Introduction 
We consider the optimal control problem: 

f T 

(la) maximize / i(x{s), u(s)) ds + <f>(x(T)) 

Jo 

over the set of trajectories (x(-), «,(•)) satisfying 

(lb) x(s) = f(x(s),u(s)), x(0)=x, x(s)eX, u{s)eU, 

for all < s < T. Here, the state space X is a subset of R™, the set of control 
values U is a subset of R m , the horizon T > and the initial condition x G X 
are given, we assume that the map u(-) is measurable, and that the map x(-) is 
absolutely continuous. We also assume that the instantaneous reward or Lagrangian 
I : X x U — > R, and the dynamics f : XxU — > R™, are sufficiently regular maps, and 
that the terminal reward <p is a map X — > RU{— oo}. The value function v associates 
to any (x,t) € X x [0,T] the supremum t) of J"' £(x(s), u(s)) ds + (f>(x(t)), under 
the constraint IjlbJ) . for < s < t. Under certain regularity assumptions, it is known 
that v is solution of the Hamilton- Jacobi equation 

Bv 3v 

(2a) = , (x,t) e X x (0,T] , 

with initial condition: 

(2b) v(x,0) = (j)(x), x e X , 

where H(x,p) = sup u6 v t(x, u) + p ■ f(x,u) is the Hamiltonian of the problem (see 
for instance |Lio82l IPS93I Lr3ar94p . The evolution semigroup 5* of J2J associates to 
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any map <f> the function v l := v(-,t), where v is the value function of the optimal 
control problem l|la|l. 

Maslov |Mas73j (see also |MS92llKM97] ) observed that the evolution semigroup 
S* is max-plus linear. Recall that the max-plus semiring, H^max? is the set IRU{ — oc}, 
equipped with the addition a © b = max(a, b) and the multiplication a <E) b = a + b. 
By max-plus linearity, we mean that for all maps /, g from X to R ma x, and for all 
A G Kmax, we have 

S t (f®g) = S t f®S t g , 
S*(A/) = AS*/ , 

where fOg denotes the map x i— » f(x)®g(x), and A/ denotes the map x i— > A<g>/(x). 
Linear operators over max-plus type semirings have been widely studied, see for 
instance IHU791 IMS921 |B(JOQ92| IKM97I l( IM0l| . 

In this paper, we introduce a new discretization method to solve the deterministic 
optimal control problem ijTJl. using the max-plus linearity of the semigroup S*. 
In I- MOD , Fleming and McEncancy introduced a max-plus based discretization 
method to solve a subclass of Hamilton- Jacobi equations (with a Lagrangian £ 
quadratic with respect to u, and a dynamics / affine with respect to it). They 
approximated the evolution semigroup S* by a max-plus linear semigroup acting 
on a finitely generated semimodule of functions. This work was pursued in [McEOll 
IMcEOOl IMcEOSbl IMcEOSa] . Another max-plus based numerical work on Hamilton- 
Jacobi equations is due to Bacaer |Bac01llBac02| . The different discretization that 
we introduce here relies on a notion of max-plus "variational formulation" , which 
originates from the notion of generalized solution of Hamilton-Jacobi equations of 
Maslov and Kolokoltsov |KM88) . )KM97I Section 3.2]. This discretization, which 
can be interpreted geometrically in terms of projections on semimodulcs, is similar 
to the classical finite element method. We shall see that the space of test functions 
must be different from the space in which the solution is represented, so that our 
discretization is indeed a max-plus analogue of the Petrov-Galerkin finite element 
method. We illustrate the method by numerical examples. We also give an error 
estimate, in dimension one, of order V At+Ax(At)~ 1 , which is the same as the order 
obtained for existing discretization methods, see |Fal87| and |BCD97I Appendix A, 
by M. Falcone] 

The present paper is only a preliminary account: the results will be detailed 
elsewhere. A first presentation of the method appeared in |Lak03j . 

2. Preliminaries on residuation and projections over semimodules 

In this section we recall some classical residuation results (see for example 
|DJLC53j . |Bir67j . |BJ72j . |BCOQ92|) , and their application to linear maps on 
idempotents semimodules fsee |LMS01I |C(iQ04| ). We also review some results of 
|CGQ96| |CGQ04| concerning projectors over semimodules. 

2.1. Residuation, semimodules, and linear maps. If (S, <) and (T, <) are 

(partially) ordered sets, we say that a map / : S — ► T is monotone if s < s' ==> 
/(s) < /( s 0- We say that / is residuated if there exists a map /" : T — » S such 
that 

f(s)<t ^ *</*(*) • 
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The map / is residuated if, and only if, for all t G T, {s G S \ f(s) < t} has a 
maximum element in S. Then, 

f\t) = max{s G S \ f(s) <t}, Vt G T . 

Moreover, in that case, we have 

(3) / o ft o / = /* and /»o/o/« = /. 

If a set /C is a monoid for a commutative idempotent law © (idempotent means that 
a ® a = a) , the natural order on /C is defined by a < 6 a(Bb ~ b. We say that 

/C is complete as a naturally ordered set if any subset of JC has a least upper bound 
for the natural order. If (JC, ©, (E>) is an idempotent semiring, i.e., a semiring whose 
addition is idempotent, we say that the semiring JC is complete if it is complete as 
a naturally ordered set, and if the left and right multiplications, L^, -R„ : JC — > /C, 
Lq(x) = ax, R^(x) = xa, are residuated. 

The max-plus semiring, R roa x 5 is & n idempotent semiring. It is not complete, but 
it can be embedded in the complete idempotent semiring M max obtained by adjoin- 
ing +oo to K niaX j with the convention that — oo is absorbing for the multiplication 
a <g> b = a + b. The map x i— » — x from K to itself yields an isomorphism from K ma x 
to the complete idempotent semiring IR min , obtained by replacing max by min and 
by exchanging the roles of +00 and —00 in the definition of M max . 

Scmimodules over semirings arc defined like modules over rings, mutatis mutan- 
dis, see |LMS01I |CGQ04| . When JC is a complete idempotent semiring, we say that 
a (right) /C-semimodulc X is complete if it is complete as a naturally ordered set, 
and if, for all u G X and A G JC, the right and left multiplications, R? : X — > X, 
v 1 ► v\ and : JC — > A", /i 1— > w/i, are residuated. In a complete semimodule <Y, 
we define, for all u, v G A", 

u\« d = = max{A G /C | uA < w} . 

We shall use semimodules of functions: when X is a set and (JC, ®, ®) is a 
complete idempotent semiring, the set of functions JC X is a complete /C-semimodule 
for the componentwise addition (u, v) 1— > m©w (defined by (tt©u)(x) = u(x)(Bv(x)), 
and the componentwise multiplication (A, it) 1— * uX (defined by (u\)(x) = u(x)®\). 

If JC is an idempotent semiring, and if X and y are /C-semimodules, we say that 
a map A : X — > ^ is additive if for all it, u G A, A(u ® u) = A(u) ® ^4(w) and 
that A is homogeneous if for all u G A and A G /C, A(uA) = A(u)A. We say that 
A is linear, or is a linear operator, if it is additive and homogeneous. Then, as in 
classical algebra, we use the notation Au instead of A(u). When A is residuated 
and v £ y, we use the notation A\v or A$v instead of A^(v). We denote by 
L(X,y) the set of linear operators from X to y. If JC is a complete idempotent 
semiring, if X, y, Z are complete /C-semimodules, and if A G L(y, Z) is residuated, 
then the map La ■ L(X,y) — > L(A,2), B h A o B, is residuated and we set 
A\C* := (Xa) j (C), for all C* G i(A,Z). 

If X and y are two sets, JC is a complete idempotent semiring, and a G /C XxY , 
we construct the linear operator A from /C r to JC X which associates to any u G /C Y 
the function Au G /C x such that ^4u(a;) = V ye y a(x,y) ® where V denotes 

the suprcmum for the natural order. We say that A is the kernel operator with 
kernel or matrix a. We shall often use the same notation A for the operator and 
the kernel. As is well known (see for instance |BCOQ92| ), the kernel operator A is 
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residuated, and 



(A\v)(y) = A A(x,y)\v(x), 



where A denotes the infimum for the natural order. In particular, when K, = M max , 
we have 

(4) (A\v)(y) = A J-A(x, y) + v(x)) = [—A* (—v)](y) , 

x£X 

where A* denotes the transposed operator JC X — » JC Y , which is associated to the 
kernel A*(y, x) = A{x, y). (In |0J, we use the convention that +00 is absorbing for 
addition.) 

2.2. Projectors on semimodules. Let V denote a complete subsemimodule of a 
complete semimodule X over a complete idempotent semiring /C, i.e., a subset of 
X that is stable by arbitrary sups and by the action of scalars. We call canonical 
projector on V the map 

(5) P V :X^X, u i-> P v (u) = max{u G V | v < u}. 

Let W denote a generating family of a complete subsemimodule V, which means 
that any element v G V can be written as v = \/{wX w \ w G W}, for some X w 6 K,. 
It is known that 

Pv(u) = V w(w\u) 

(see for instance |CGQ04| ). If B : U — > X is a residuated linear operator, then the 
image im B of B is a complete subsemimodule of X, and 

(6) P imB = BoBK 

The max-plus finite element methods relies on the notion of projection on an im- 
age, parallel to a kernel, which was introduced by Cohen, the second author, and 
Quadrat, in |CGQ96| . The following theorem, of which Proposition [3 below is an 
immediate corollary, is a variation on the results of |CGQ96[ Section 6] . 

Theorem 1 (Projection on an image parallel to a kernel). Let B : U — * X and 

C : X — > y be two residuated linear operators. Let LI§ = B o (C o B)^ o C . We 

have Llg = LI B o n c , where LI B = B o B^ and U c = o C. Moreover, Llg is a 
projector ((LI§) 2 = LT§), and for all x G X: 

H%(x) = max{y e imB \ Cy < Cx}. 

The results of |CGQ96| characterize the existence and uniqueness, for all x G X , 
of y £ imB such that Cy = Cx. In that case, y = LI^(a;). 

When JC = IR m ax, and C : M max — > K max is a kernel operator, U c = C'oC has 
an interpretation similar to ©: 

n C («) = C« o C(«) = -F imC .(-^) = P-imO*(«) , 
Y 

where — imC* is thought of as a R ro i n -subsemimodule of M min , so that, 

P-imC*(v) = mm{w G — imC* | w > v} . 

x 

where < denotes here the usual order on M , since the natural order of M m j n is the 

u x 

reverse of the usual order. When B : M max — > R max IS a l so a kernel operator, we 
have 

II B = PimB P-imC* ■ 
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This factorization will be instrumental in the geometrical interpretation of the finite 
element algorithm, see Example 1101 below. 

3. The max-plus finite element method 

3.1. Max-plus variational formulation. We now describe the max-plus finite 
element method to solve the optimal control problem l|la(l . Let S* and v l be defined 
as in the introduction. Using the semigroup property S t+t = S* o 3* , for t, t' > 0, 
we have the recursive equation: 

(7) v t+At = S* A V, t = 0,Af,--- ,T- At 

with v° = 4> and At = ^ , for some positive integer TV. Let W be a R max -semimodule 
of functions from X to R max such that <j> £ W and for all v £ W, t > O^v £ W. 
We suppose given a "dual" semimodule Z of "test functions" from X to M m ax- The 
max-plus scalar product is defined by (u \ v) = sup^g^ u(x) +v(x), for all functions 
u, v : X — > R, with the convention that — oo is absorbing for the addition +. We 
replace {7} by: 

(8) (z | v t+At ) = (z | S ,At w t ), Vz£Z, t = 0,At,...,r- At , 

with v At , . . . ,v T £ W. Equation © can be seen as the analogue of a variational 
orweak formulation. Kolokoltsov and Maslov used this formulation in |KM88| and 
|KM971 section 3.2] to define a notion of generalized solution of Hamilton- Jacobi 
equations. 

3.2. Ideal max-plus finite element method. We consider a semimodule Wh C 
W with generating family {wi}i<i< p . We call finite elements the functions Wi. We 
approximate v t by £ Wh, that is: 

«* - v h= V w t A* , 

where A* £ M ma x- We also consider a semimodule Z/j £ Z with generating family 
{■Zj}i<j<<j- The functions z\, ■ ■ ■ , z q will act as test functions. We replace JSJ by 

(9) ( z j I V h +At ) = ( z i I sAtv l)> Vl<i<g , 

for t = 0, At, • ■ • , T — At, with = h ~ and £ W h , t = 0, At, • ■ • , T. 

Since Equation (j9j) need not have a solution, we look for the maximal subsolution, 
i.e. the maximal solution v h +At £ Wh of 

(10a) feK +At ) < fe|5 A X> Vl<i<g . 

We also take for the approximate value function at time the maximal solution 
vl £ W h of 

(10b) v° h < v° . 

Let us denote by Wh the max-plus linear operator from Mn iax to W with matrix 
Wh = co\(wi)i<i< p , and by the max-plus linear operator from W to M' nax whose 
transposed matrix is Zh = col(^j)i<j< 9 . This means that WhX = V \<i< p Wi\i for 
all A = (Ai)i=i,..., p £ Mf nax , and (Z^v)j = (zj \ v) for all v £ W and j = 1, . . . , q. 
Applying Theorem ^ to B = Wh and C = Z^ and using WV = imWh, we get: 
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Proposition 2. The maximal solution v k +At 6 Wh of IjlOafl is given by v k +At = 
S At v k , where 

S At = Il§ h o S At . 

Proposition 3. Let v h £ Wh be the maximal solution of 1)10(1 . for t = 0, At, . . . , T . 

Then, for every t = 0, At, ...,T, there exists A* € I^max suc h that v\ = Wf l X t . 
Moreover, the maximal A* satisfying these conditions verifies the recursive equation 

(Ha) A t+At = (ZtWhUZtS^W^) , 

with the initial condition: 

A = W h \</> . 

Proof. Since v\ <E Wh, vt = WhX f , and the maximal A* satisfying this condition is 
A* = Wl{v h ), for all t = 0, At, . . . , T. Since w° is the maximal solution of HlObJI . then 

by © and ©, v° h = P Wh {4>) = W h oWl{4>), hence A = Wj|oW),oW*(0) = ^ (</>)■ 
Let t = 0, . . . , Ta^. Using Proposition |2 Theorem and the property that 

(/ ° .9)" = ff" /" f° r a ll residuated maps / and g, we get 

A t+At = wloU^oS^iWhX 1 ) 

= W{ o W h o Wjg o {Z*J o Z* o s*\w h \*) 

= Wi o (^)« o Zl o 5 At (M/ /l A t ) 

= (ZlWh^iZtS^WhX 1 ) . 
which yields Hlla(l . □ 

The maps A h := ZjJW}, : Mf llax -> and B ft := Z*S At W^ : -> M? nax 

are max-plus linear operators, and the entries of their corresponding matrices are 
given, for 1 < i < p and 1 < j < q, by: 

(12) {A h )ji = {z 3 I Wi) 

(13) {B h ) 3l = ( Zj I S At Wl ) 

(14) = ((S*) At z, I wi) , 

where 5* is the transposed semigroup of S, which is the evolution semigroup associ- 
ated to the optimal control problem in which the sign of the dynamics is changed. 
The ideal max-plus finite element method can be summarized as follows: 

(1) Choose At = jj and the finite elements (w>i)i<i< p and (Zj)i<j<q, 

(2) Compute the matrix Ah by 1|12|) and the matrix P>h by (|13|) or by (|14fl . 

(3) Compute A = W h \(f> and v° h = W h X°. 

(4) For t = At, 2At, ...,T, compute A* = A /l \(S h A t " At ) and v\ = W h A*. 

Then, v k approximates the value function at time t, v . 

The recursion A f = A/ l \(_B^A t_At ) may be written explicitly as 

A* = mm ( - (A*)* + max ((B h ) jk + A^ A *)) , for 1 < i < p . 

Observe that this recursion may be interpreted as the dynamic programming equa- 
tion of a deterministic zero-sum two players game, with finite action and state 
spaces. 
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In order to implement this method, we must specify how to compute the entries 
of At and Bh in l|12|) and (|13|l or Ijl4jl . In some cases, these computations can be 
done analytically. Computing Ah from l|12|) is an optimization problem which may 
be solved by standard algorithms. We shall discuss in the following section the 
approximation of Bh- 

3.3. Effective max-plus finite element method. We first discuss the approx- 
imation of S At w for every finite clement w. The Hamilton- Jacobi equation l|2a(l 
suggests to approximate S At w by the function [S At w]~ such that 

(15) [S At W r(x) = w(x)+AtH(x, forallzeX. 

ax 

Let [S At Wh}~ denotes the max-plus linear operator from K^ lax to W with matrix 
[5 A *^]~ = col^S^WtY )i< l < P , which means that 

[S M W h )~\ = V [S M Wi]~^ 

l<i<p 

for all A = (Xi)i<i< p £ KJ| 1M . The above approximation of S At w yields an approx- 
imation of the matrix Bh by the matrix B^ := Z^[S At Wh]~ , whose entries are 
given, for 1 < i < p and 1 < j < q, by: 

(Bh )ji = sup (zj{x) + Wi(x) + AtH(x, -7T 1 )) ■ 

x£X OX 

Thus, computing B^ requires to solve an optimization problem, which is nothing 
but a perturbation of the optimization problem associated to the computation of 
Ah- We may exploit this observation by replacing B^ by the matrix B^~ with 
entries 

(16) (BDii = (zj\wi)+At sup H{x,-p-) , 

for 1 < i < p and 1 < j < q. Here, argmaxjzj + wj denotes the set of x such 
that Zj(x) + Wi(x) = (zj | Wi). When this set has only one element, (|16() yields a 
convenient approximation of Bh- 

Of course, Wi must be differentiable for the approximation l|15|l to make sense. 
When Wi is non-diffcrentiable, but Zj is differentiable, we may approximate (Bh)ji 

by 

dz ■ 

sup(zj(x)+AtH(x,-—^)+w i (x)) , 

x£X OX 

using the dual formula l|14f) . We may also use the dual formula of i|16[l . where 
is replaced by — 

3.4. Comparison with the method of Fleming and McEneaney. Fleming 
and McEneaney proposed a max-plus based method |FM00| . which also uses a 
space Wh generated by finite elements, Wi, . . . , w p , together with the linear formu- 
lation (J7J). Their method approaches the value function at time t, v f , by Wh^*, 
where Wh = co\(wi)i<i< p as above, and ^* is defined inductively by 

(17a) M ° = W h \(f> 

(17b) M t+At = (W h \{S At W h ))^ , 
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for t = 0, At, . . . , T — At. This can be compared with the limit case of our finite 
element method, in which the space of test functions Zh generates the set of all 
functions. This limit case corresponds to replacing by the identity operator 
in (|lla[) . so that 

(18) A* +At = W h \{S M W h \*) • 

Proposition 4. Let (//*) be the sequence of vectors defined by the algorithm of 
Fleming and McEneaney, (|17(l : let (A*) be the sequence of vectors defined by the 
max-plus finite element method, in the limit case (|18|) ; and let v denote the value 
function at time t. Then, 

W h rf < WhX 1 < v* , for t = 0, At, . . . , T . 

Sketch of proof . This can be proved by induction, by using the residuation inequal- 
ity WlS At W h X > (W h \(S At W h ))X, which holds for all vectors A, together with 
the monotonicity of the operators arising in the construction of A* and fi . □ 

An approximation of (|17b|) using formulae of the same type as l|15f) is also dis- 
cussed in |MH99| . An experimental comparison will appear elsewhere. 

4. Error analysis 

The following general lemma shows that the error of the finite element method 
is controlled by the projection errors, HIL^- u* — u*||oo and ||n z ' i t; t — u*||oo) an d by 
the approximation errors, || [£ At u>i]~ — S At Wi\\oo, and \(B^~)ji — (B^)ji\- 

Lemma 5. For t = 0, At, ■ ■ ■ ,T, let v l be the value function at time t, and v h be its 
approximation given by the effective max-plus finite element method, implemented 
with the approximation of Bh, given by (|16|) . We have 

\\vl-v T \\ 00 <(i + ^-)( sup (iin^^-^iioo + nn^y-^iu) 

At \ t =0,At,...,T 

+ max || [S A V]~ " ^ At ^Hoo + max \{B~ h ~ )ji ~ )*|) . 
i<i<p ^SjSq ' 

l<i<p 

The proof of this lemma uses the fact that projectors over max-plus semimodules 
are non-expansive in the sup-norm. 

To state an error estimate, we make the following assumptions: 

- (HI) The semigroup preserves the set of i-semiconvex functions, for some 
c> 0. 

- (H2) f : X x U — > W l is bounded and Lipschitz continuous with respect to 
x: 

3L f >0, Vx,yeX, \f(x,u)-f(y,u)\<L f \x-y\ Vu£U, 
3M/>0, Vx,yeX, \f{x,u)\<M f . 

- (H3) £ : X x U — > M is bounded and Lipschitz continuous with respect to 
x: 

f \£(x,u)~e{y,u)\<Li\x-y\ Vx,y eX,ueU, 
\ \£(x,u)\<M h Mx,yeX,ueU. 

- (HA) (f> : X — *■ K is bounded and Lipschitz continuous: 

\(f>(x) -<f>(y)\ < L^\x-y\ Vx,yeX. 
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Recall that a function / is ^-semiconvex if f(x) + ^x 2 is convex. Spaces of semi- 
convex functions were already used by Fleming and McEneaney |FM00) . 
We shall use the following finite elements. 

Definition 6 (Lipschitz finite elements). Assume that X is an interval of R. We 
call Lipschitz finite element centered at point x € X, with constant A > 0, the 
function w(x) = — A\x — x\. 

Definition 7 (Quadratic finite elements). Assume that X is an interval of R. We 
call quadratic finite element centered at point x G X, with Hessian - > 0, the 
function w(x) = — ^{x — x) 2 . 

The family of Lipschitz continuous finite elements of constant A generates, in 
the max-plus sense, the semimodule of Lipschitz continuous functions of Lipschitz 
constant A. When X = R, the family of quadratic finite elements with Hessian i 
generates, in the max-plus sense, the semimodule of lower-semicontinuous i-semi- 
convex functions. 

Theorem 8. Let X = [—b,b] C R. We make assumptions (H1)-(H4), and assume 
that there exist L > such that the value function at time t, v , is L-Lipschitz 
continuous and ^-semiconvex for all t > 0, with the same constant c as in (HI). 
Let us choose quadratic finite elements w; of Hessian -, centered at the points of 
the regular grid (ZAx)n [— (b + cL), (b + cL)]. Let us choose, as test functions Zj, the 
Lipschitz finite elements with constant A> L, centered at the points of the regular 
grid (ZAx)n [—6, b] . For t — 0, At, . . . , T , let v l h be the approximation of v l given by 
the effective max-plus finite element method, implemented with the approximation 
Bfr ~ of Bh- Then, there exists a constant K > such that, for At small enough, 

H-v T \\ 00 < I <(VAt+^) . 

A variant of this theorem, with a stronger assumption, is proved in |Lak03| . We 
shall give elsewhere the proof of Theorem [SI 

5. Numerical results 

Example 9 (Linear Quadratic Problem). We consider the case where U = R, 
X = R, 

i{x,u) = -(^| 2 + ^-), f(x,u)=u, and = . 

We obtain H(x,p) — — f |x| 2 + \- We choose quadratic finite elements Wi and 
Zj of Hessian 1, centered at the points of the regular grid (ZAx) n [— L,L]. We 
represent in Figure ^ the solution given by our algorithm in the case where T = 5, 
At = Ax = 0.05, a = 0.3 and L = 10. The computations were performed using the 
max- plus toolbox of Scilab |Plu98| . 

Example 10 (Distance problem). We consider the case where T = 1, <f> = 0, 
X = [-l, 1], Z7 = [-1,1], 

<(*,«) = ( _1 if xe( r M) ; and f(x, u) = h if ^(-m). 

V ' [0 if are {-1,1}, (0 if ie{-l,l}. 

Consider first quadratic finite elements u>i and Zj of Hessian 1, centered at the 
points of the regular grid (ZAx) n [—1,1]. In Figure [3 we represent the solution 
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Exact solution 



Approximated solution 



10 



Figure 1. Max-plus approximation of a linear quadratic control 
problem (Example 



given by our algorithm in the case where At = 0.05, Ax = 0.0125 and c = 1.2. Since 
U Zh is a projector on a subsemimodule of the M min -scmimodulc of — i-scmiconcavc 
functions, and since the solution is not — ^-semiconcave for any c, the error of 
projection ||n z,i (u*) — u'||oo does not converge to zero when Ax goes to zero, which 
explains the magnitude of the error. 




— Exact solution 
Approximated solution 



Figure 2. A bad choice of test functions for the distance problem (Example 1101 

To solve this problem, it suffices to replace the test functions Zj by the Lipschitz 
finite elements with constant A > 1, centered at the points of the regular grid 
(ZAx) n [—1,1]. This is illustrated in Figure [3] in the case where At = 0.05, 
Ax = 0.0125, c = 1.2 and A = 1.1. 

The next two examples are inspired by those proposed by M. Falcone in |BCD97| . 

Example 11. We consider the case where T = 1, $ = 0, X = [—1, 1], U = [0, 1], 
£(x,u) = x and f(x,u) = —xu. The optimal choice is to take u* = whenever 
x > and to move on the right with maximum speed (it* = 1) whenever x < 0. 
For all t £ [0,T], the value function is: 

, (xt if x > 

v(x,t) = < 

I — e ) otherwise. 
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Exact solution 

Approximated solution 



Figure 3. A good choice of test functions for the distance problem ( Example 110(1 



We choose quadratic finite elements u>; of Hessian - and Lipschitz finite elements Zj 
with constant A > 1. We represent in Figure 0] the solution given by our algorithm 
in the case where T = 1, At = 0.05, Aa: = 0.02, A = 1.3 and c = 1.4. 




0.1 - 



Exact solution 
Approximated solution 



Figure 4. Value function and its max-plus approximation (Example II If) 



Example 12. We consider the case where T = 1, $ = 0, X = [— 1, 1], U = [-1, 1], 
£(x, it) = —3(1— \x\) and f(x, u) = u(l— \x\). The optimal choice is to take u* = — 1 
whenever x > and it* = 1 whenever x < 0. For all i e [0, T], the value function 
is: 

v(x,t) = -3(1 - |x|)(l - e _t ). 

We choose quadratic finite elements u>^ of Hessian i and Lipschitz finite elements 
Zj with constant ^4. We represent in Figure the solution given by our algorithm 
in the case where T = 1, At = 0.05, Ax = 0.02, A = 2 and c = 1.1. 
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Exact solution 
Approximated solution 



-2.0 



-1.0 
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Figure 5. Value function and its max-plus approximation (Example El 
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